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' Recent numerical evidence suggests that a mass spectrum of primordial 

J> ■ black holes (PBHs) is produced as a consequence of near critical gravitational 

collapse. Assuming that these holes formed from the initial density perturba- 
O , tions seeded by inflation, we calculate model independent upper bounds on the 

mass variance at the reheating temperature by requiring the mass density not 

' exceed the critical density and the photon emission not exceed current diffuse 

gamma-ray measurements. We then translate these results into bounds on 

• the spectral index n by utilizing the COBE data to normalize the mass vari- 

ance at large scales, assuming a constant power law, then scaling this result 
to the reheating temperature. We find that our bounds on n differ substan- 
tially (5n > 0.05) from those calculated using initial mass functions derived 
under the assumption that the black hole mass is proportional to the horizon 
mass at the collapse epoch. We also find a change in the shape of the diffuse 
gamma-ray spectrum which results from the Hawking radiation. Finally, we 
study the impact of a nonzero cosmological constant and find that the bounds 
on n are strengthened considerably if the universe is indeed vacuum-energy 
dominated today. 
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I. INTRODUCTION 



Primordial black holes (PBHs) arise naturally in most cosmologies. Perhaps the least 
speculative mechanism for PBH formation comes from the collapse of overdense regions in 
the primordial density fluctuations that gave rise to structure in the universe ffl. Thus 
PBHs carry information of an epoch about which we know comparatively little, and are 
a very useful tool for restricting theories of the very early universe, especially within the 
context of an inflationary scenario. In this paradigm the density spectrum is a consequence 
of the quantum fluctuations of the inflaton field, and can in principle be calculated given 
an underlying model. Thus, the non-observation of the by-products or of the effects of the 
energy density of these PBHs constrains the underlying microscopic theory. 

The simplest bound that can be extracted from PBH formation is generated by insisting 
that VtpBH < 1- Other bounds may be derived by studying the consequences of their 
evaporation. Given that black holes evaporate at a rate proportional to their inverse mass 
squared 0, the phenomenological relevance of a PBH will depend upon its initial mass. 
Smaller mass PBHs (10 9 < < 10 13 g) will alter the heavy elements abundances || as 
well a distort the microwave background whereas PBHs with larger masses will affect 
the diffuse gamma-ray background 

The net evaporation spectrum from a collection of PBHs will depend on the initial mass 
distribution, which in turn depends upon the probability distribution for the density fluc- 
tuations. Prior to the emergence of the COBE data, bounds from PBHs were calculated 
assuming a Harrison-Zel'dovich spectrum \Sk\ 2 oc k, with an unknown normalization. This 
lead to the famous "Page-Hawking" bound and all its subsequent improvements [|5]|| . How- 
ever, assuming that the distribution is Gaussian, we can now relate the mass variance at 
the time of formation with the mass variance at large scales today, if we know the (model- 
dependent) power spectrum. None the less, bounds on n have been derived within the class 
of models where the power spectrum is given by a power law \Sk\ 2 oc k n over the scales of 



interest [p|-|T2|j. Indeed, in this way it follows that for the scale-invariant Harrison-Zel'dovich 
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spectrum PBHs have too small a number density to be of any astrophysical significance. 
However, observations (see references in |J) now seem to favor a tilted blue spectrum with 
n > 1 (in CDM models) with more power at smaller scales, and thus the bounds derived 
from black hole evaporation can be used to constrain the tilt of the spectrum. 

In this paper we revisit the aforementioned bounds in light of recent calculations which 
indicate that a spectrum of primordial black hole masses are produced through near critical 
gravitational collapse. As was pointed out by Jedamzik and Niemeyer |13[ , if PBH formation 
is a result of a critical phenomena, then the initial mass function will be quite different then 
what was expected from the classic calculation of Carr [Q. In particular, the PBH mass 
formed at a given epoch is no longer necessarily proportional to the horizon mass. The 
resulting difference in the initial mass function leads to new bounds, which is the main 
thrust of this paper. In particular, we revisit the density bounds and the bounds derived 
from the diffuse gamma-ray observations. We will first derive bounds on n in the class of 
models where the power spectrum is a power law over the scales of interest. We then relax 
this assumption and instead place bounds on the mass variance at the formation epochfj 
These bounds can then be applied to a chosen model by extrapolating the variance today 
to the formation epoch according to the appropriate power spectrum. 

Before continuing to the body of this work, we would like to point out that primordial 
black holes have also played a large role in attempting to explain various data. PBHs can 
serve as significant cosmological flux sources for all particle species via Hawking radiation || . 
Thus it is very tempting to postulate that present-day observed particle fluxes of unknown 
origin are a consequence of PBH evaporation. However, to predict these fluxes, or model 
them realistically, we need to know the mass distribution of black holes, since their emission 
spectra are determined by their temperature (or inverse mass) . If we assume that the black 



x As will be discussed later, PBH formation is dominated by the earliest formation epoch if they 
form via critical collapse. 
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holes of interest were formed from initial density inhomogeneities generated in an inflationary 
scenario (which is usually assumed), then the black holes are either tremendously over- 
abundant or completely negligible. To get a phenomenologically interesting quantity of 
PBHs thus requires an extreme fine-tuning, as will be demonstrated below. Succinctly, this 
fine-tuning arises because the PBH number density is an extremely rapidly varying function 
of the spectral index n. Thus, without even analyzing the details of the spectral profile, 
explaining unknown fluxes via PBH evaporation is far from compelling. 



II. THE INITIAL MASS FUNCTION 



Carr |14| first calculated the PBH spectrum resulting from a scale-invariant Harrison- 
Zel'dovich spectrum up to an overall normalization. Subsequently Page and Hawking cal- 
culated a bound on the normalization by calculating the expected diffuse gamma-ray spec- 
trum from these PBHs j5]|| . However, using the COBE measurements of the temperature 
anisotropies translated into density fluctuations (within a CDM), the overall normalization 
can be determined. For a scale-invariant spectrum, no significant number density of PBHs 
is generated. However, for a tilted blue power spectrum with more power on smaller scales, 
a larger number density of PBHs is expected. 

Given an initially overdense region with density contrast Si and radius R at time ti (using 
the usual comoving coordinates), analytic arguments predict |],[14| a black hole will form if 

1/3 < St < 1 . (1) 

The lower bound on the density contrast comes from insisting that the size of the region at 
the time of collapse be greater than the Jeans length, while the upper bounds come from 
the consistency of the initial data with the assumption of a connected topology. 
The first calculations of the PBH mass distribution assumed the relation 

M hh ~ 7 3 / 2 i\4, (2) 
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where 7 determines the equation of state p = 7p and Mh is the horizon mass when the 
scale of interest crossed the horizon. Recently, numerical evidence suggests that near the 
threshold of black hole formation, gravitational collapse behaves as a critical phenomena 
with scaling and self-similarity A scaling relation of the following form was found 

M bh (5) = kM h (5-5 c y, (3) 

where p is a universal scaling exponent which is independent of the initial shape of the 
density fluctuation. It was later shown [13] that such scaling should be relevant for PBH 



formation. Indeed, in Ref. [[HJ the authors found relation (|3]) to hold for PBH formation 
when the initial conditions are adjusted to be nearly critical. They found the exponent to 
be p ~ 0.37. They also found that for several different initial density shapes, 5 C ~ 0.7, which 
is significantly larger than the analytic prediction of 1/3 found by requiring that the initial 
overdensity be larger than the Jeans mass. 

Given Eq. ([3]), calculating the initial PBH mass distribution becomes analytically cum- 
bersome, since in principle one needs to sum over all epochs of PBH formation. However, 
we would expect that the initial mass function would be dominated by the earliest epoch 
of formation if we assume a Gaussian distribution with a blue power spectrum, since for 
larger scales the formation probability should be suppressed. This expectation was tested 



in Ref. ||17|1 , where the authors used the excursion set formalism |L8|] to calculate the initial 
mass function allowing for PBH formation at all epochs. The authors found that it was 
approximately true the that earliest epoch dominates, for the conditions of interest to us. 



This simplification allows us to derive quite easily the initial mass distribution [fL3|| . We 
assume Gaussian fluctuations (the effects of non-Gaussianity will be briefly discussed in 
Sec. 0) and define the usual smoothed density contrast 

6 R (x) = Jd 3 y6(x + y)W R (y) , (4) 

where 8(x) = (p(x) — Pb)/Pb, and pt, is background energy density. Wr is the window 
function with support in a region of size R. The probability that a region of size R has 
density contrast between S and S + dS is given by 



P{R,5)d5 



exp 



8 2 



d5 



(5) 



where <jr is the mass variance for a region of size R, a 2 R = {S R (x))/R 3 . 

Then using Eq. (^), the physical number density of PBHs within the horizon per loga- 
rithmic mass interval, at the formation epoch, can be written as 



d log M bh 



dlogMi 



bh 



/ Z7T a p 



\kM H ) 



VP 



exp 



2a 2 



S c + 



( M x 



bh 



\kM 



H 



(6) 



where Vh is the horizon volume the the epoch of PBH formation. We assume prompt reheat- 
ing, and therefore take the formation epoch to be the time of reheating, which corresponds 



to the minimum horizon mass 17 



Assuming that the power law spectrum holds down to the scales of the reheat temper- 
ature, then <7jj oc R~( n+3 \ We can then relate the mass variance today ctq to the mass 



variance at the epoch of PBH formation a(M R ), using [10,17 



a 2 (M H ) = o\ 



M eq \^~ n ^ (M H 



Mn 



(l-n)/2 



M, 



(7) 



where 



M, 



H 



M ' Teq 



2 / T \ 3/2 



(8) 



-TrhJ \T eq ) 

M is the mass inside the horizon today, Trh is the reheat temperature, and T eq is the 
temperature at radiation-matter equality. This relation is essential to connect present-day 
density fluctuations to those of much earlier times. Two important assumptions underly 
this useful relation: First, n is taken to be constant over all the scales of interest. Second, 
the universe was assumed to be radiation dominated until the temperature dropped below 
T eq ~ 5 eV (matter-radiation equality), and then matter dominated thereafter. 



From the COBE anisotropy data, the mass variance can be calculated [19,10 



ct = 9.5 x 10 



-5 



(9) 
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Using this result we can then calculate the physical number density per unit mass interval 
at T = T RH 



dn, 



bh 



VT 1 



dM, 



bh 



2ira(M H )M bh p 



y l l p exp 



s c + y 



2a 2 (M j 



ii) 



(10) 



where 



bh 



M bh T RH 



kM H 0.301A;^ 1/2 M| 



(11) 



in 



and Mpi is the Planck mass. The physical number density at time t is simply Eq. ( |T0D 
rescaled by a ratio of scale factors, that can be written as 



dn bh . _ ( T(t) \ dn bh 



(12) 



dM bh y \Tr„J dM bh 

Finally, we should note that relation (|3|) is only valid for 5 ~ 8 C . Thus we expect that we 
may integrate over 5 with small errors, as long as the width of the Gaussian is sufficiently 
small. In particular, our results should be trustworthy provided a ^ 1, which implies n 
should not exceed the maximum 

21og(a 2 ) 



n r , 



1 + 



(13) 



\og(T eq T Q /Tl H ) ■ 

We will see that n does not exceed this maximum value for all of the bounds we consider 



III. BOUNDS FROM O < 1 

Let us now calculate the total energy density in PBHs. We assume the "standard cosmol- 
ogy" where the universe began in an inflationary phase, reheated, was radiation dominated 
from the reheating period until matter-radiation equality, and then has been matter dom- 
inated. The contribution of a PBH with a given initial mass, M b h, to the energy density 
today will depend upon its lifetime. The time-dependent PBH mass M(t) is given by [§] 

1/3 

(14) 



M{t) = M„ 



t 
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where M* is the initial mass of a PBH which would be decaying today, ~ 5 x 10 14 g. 
It is a good approximation to assume that the black hole decays instantaneously at a fixed 
decay time, td, which we use in the following. 

There are two components to the PBH density bounds that we can calculate. The 
first is the total energy density of the PBHs that have not decayed by a given time t. 
The second is the total energy density of the products of PBH evaporation. The sum 
of these components must be less than the critical density Q p bh,eva P + Q P bh < 1, at any 
time. The evaporated products of PBHs, in particular photons, could break up elements 
during nucleosynthesis, disrupting the well-measured elemental abundances. This and other 
processes during nucleosynthesis provide additional bounds on the density of PBHs || that 
we do not discuss here. 

The simplest bound comes from the density of PBHs that have not decayed by time t, 



where T(t) is the temperature of the universe at time t, and M±(t) is the initial PBH mass 
that has just completely evaporated by time t, 

/t\ 1/3 , , 



J Ml 2 (£) is a dimensionless weighted integral over the PBH mass spectrum between Mi to 
M 2 , normalized to the total density ptotitRii), 



where dn^h/dM^ is given by Eq. fllQ|). We can use the above to trivially compute the ratio 
of the PBH density to the critical density, Q(t). In particular, we need only compute the 
density ratio at three relevant epochs: immediately after reheating t = t R H, at matter- 
radiation equality t = t eq , and present-day t = t . The density ratios areQ 

2 Note that £l(t RH ) = I^ max {0) is often denoted by (3{t RH ). 




(15) 




(16) 




(17) 
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ft(fflff) = k ) max (0) (18) 

«W = ^aSt(o) (19) 

eg 

fi(to) = ^/jSr (0) • (20) 

eg 

where M max corresponds to 5 = 1, which is of order the horizon mass at reheating. The 
integral should be independent of the upper limit M max if we are to trust our results. 

Making the conservative approximation that all the PBH decay products are relativistic, 
the contribution to the density ratio of the products of PBH evaporation that has occurred 
up until today can be written as 

fW*(*o) = ^ f^Y V*(3/2) + ^t( 2 ) • ( 21 ) 

eg \ eg / -'eg 

In Fig. 0, we show the upper limit on n as a function of Trh coming from bounding 
Qpbh,evap(to) + Qpbh{to) < 1 (solid line). For larger values of the reheat temperature we 
get a more stringent bound by imposing the constraint flpbh{tRn) < 1, simply because as 
Trh is increased more of the black holes will have decayed at an earlier epoch. Given 
that most of the energy of the decay products resides in radiation, the effect on ft p bh{to) is 
diminished due to the redshifting. This new bound is given by the dotted line in Fig. [I]. 

If we assume that the PBH leaves behind a Planck mass remnant, then we have additional 
bounds which become important for very large reheat temperature [p0~|,2l|, |9| . The best bound 
in this case comes from calculating fl r emnant(t e q) which is given by 

T^RH ^^Pl A/ 

^remnantiteq) = ~M~^® ^ ' (22) 

The bound in this case is shown as the dashed line in Fig. [I], and is the best bound at the 
largest values of the reheat temperature. 

IV. BOUNDS FROM DIFFUSE GAMMA-RAYS 

For a certain range of Trh we can improve our bounds on n from diffuse gamma-ray 
constraints. The present day flux is determined by convoluting the initial mass function 
with the black hole emission spectrum 



27rexp(87rxf-(-l) 2s ' (23) 
where s is the spin of the emitted particle, x = cu (t)M(t)/M|„ u(t) is the frequency and 
M(t) is the PBH mass at the time t of emission. T s (x) is the absorption coefficient and 
may be written as [u)(t)] 2 a s /iT. a s is the absorption cross section and is calculated using the 
principle of detailed balance. 

The values for a s were calculated some time ago by Page [p2,p3]. Let us consider how 



a s behaves for massless particles. At large values of x, a s performs small oscillations about 
the geometric optics limit of a g = 27irM 2 / 'Mp X . As x approaches zero, a s goes to zero for 
s = 1/2, 1 but goes to a constant value for s = 0. We will use the approximation 

T s (x) = (56.7, 20.4)x 2 /tt for s = (|, 1). (24) 

This approximation is poor at low energies, as it is in error by 50% at x — 0.05. However, as 
we shall see, the contribution to the spectrum of interest is greatly peaked at x ~ 0.2. The 
case of strongly interacting particles is complicated by the hadronization process. There is 
a large contribution coming from pion decay, however, given the extreme sensitivity of the 
flux to the value n, the effect on the bound is negligible. 

It has been recently suggested that the self-interactions of the emitted particles 

will induce a photosphere, thus distorting the spectrum considerably from Eq. fl23f ). It was 
suggested that two types of photospheres should form. A QCD photosphere^ generated by 
parton-parton interactions as well as a QED photosphere generated by electron-positron- 
photon interactions. This is idea has been tested more quantitatively via a numerical solution 
of the Boltzmann equation [p6fl . Again, while this effect may change the spectrum, especially 



at higher energies, it is irrelevant as far as the bound on the spectral index is concerned. 
The flux measured today is given by 

dJ 1 r^o r dribh 



dun 



Lj\ t{1 + z) j dMih ^_ {t)m , (25) 



3 In the case of QCD what is meant by "photosphere" is a quark-gluon cloud. 
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where dribh/dMbh is evaluated at time t using Eq. (|TJ), to is the age of the universe, U is the 
time of last scatter, and f(x) is the instantaneous emission spectrum given above with 



uj{t)M{t) = u (l + z) M ^ 



M 2 Pl 



M 2 Pl 



Mi 



hh 



t 



1/3 



(26) 



J t . 

The integral over t is cut off at early times, since at redshifts above z = z ~ 700 the optical 
depth will be larger than unity due to either pair production off of matter or ionized matter 
p7| . Those processes will degrade the energy below the window we are interested in. 
This integral may be rewritten in the more illuminating form 



dz 



dJ 1 M% x r*o 
H7 ~ 4^ (u; M*) 3 Jo H (l + zfl 2 ju 



dxx 2 a 2 f(x) 



dn bh (x,z) 



dM, 



bh 



(27) 



where 



a 



M(t) 



-3/2 



X 3 M% 



1/3 



(28) 



Mi, \ y ~ ' ~> [(l + ^VoM,] 3 J 

Let us study the qualitative behavior of the above integral as a function of uo at fixed 
n and Trh- The x integration is controlled by the Boltzmann factor in f(x). Indeed, a 
little manipulation shows the the integrand is highly peaked near x ~ 0.2. Furthermore, 
the u dependence in a~ 2 is almost completely canceled by the u dependence in the factor 
M h ~^y x l p oc ct ( - 1 / p ~ 1 ) in dribh/dMbh- Thus the uj dependent part of the integrand may be 
written as 



dJ 

dujQ 



oc uj exp 



8 C + aTt/ga 1 /' 1 



2a 2 (M 



H 



(29) 



where a p = M+gV 2 /(0.301A; Mpj), and the only uj dependence in the exponential is through 
a. If for now we assume that the dominant contribution the higher energy photons comes 
from recent decays (z ~ 0), and most of the support for the x integral comes with x ~ 0.2, 
then a simplifies to 

1/3 



a 



1 + 



0.2M& 

. woM* , 



(30) 
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As uj gets larger than 0.2Mj, ; /M* ~ 100 MeV, a becomes independent of luq and therefore 
the flux behaves as dJ/dujQ oc cu^ 3 . For lower energies we can make the approximation 
a ~ 0. 2Mpi/(woM ic )j and we would expect that at some point the u dependence in the 
exponential will begin to dominate such that the flux should begin to rapidly decrease as we 
go to lower photon energies. The energy at which the flux turns over is determined by the 
competition between the two terms 5 C and aT^a 1 ^ in the exponential, which is set by the 
reheat temperature. As we lower the reheat temperature the position of the kink moves to 
lower energies. If the reheat temperature is higher than Trh ~ 10 9 GeV, however, the peak 
will stay around 100 MeV, since at these temperatures the second term in the exponential 
will always dominate. Indeed, we expect the position of the fall off to be near 

, klnk . mi n(l00MeV, ; 30 f;g^ ). (31) 

Figure |2] shows the flux for fixed n for a few different reheat temperatures. The position of 
the kink is well tracked by Eq. (|31]). Note however that the flux does not fall off exponentially 
at energies below the kink. This is because as we go to lower energies we pick up more of a 
contribution from higher redshifts. 

It is interesting to contrast this behavior with the flux calculated assuming that the mass 
of a PBH formed at a given epoch is proportional to the horizon mass at the time of collapse. 



In Refs. flliyi2|l the authors calculated an initial mass function following the Press-Schecter 
formalism, summing over all epochs and assuming the relation M^h — 7 3 ' 2 M# at each epoch. 
They found 

dn bh n + 3 [2 7/4 -1/3 .-5/2 -1 / 7 2 \ / QoA 
where pi and are the energy density and horizon mass at Tr# and 

( M - V" n)/4 ™ 

aH = ao {^Hfo • (33) 



This result reduces to the initial mass function first computed by Page ]TJ| for the Harrison- 
Zel'dovich spectrum with n = 1 and dnbh/dMbh oc M h ^ 2 . The u dependence of this result 
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arises only through M(t) given by Eq. (|28|). Using this initial mass distribution in Eq. 
we expect, as in the previous case, dJ/dco oc u^ 3 for larger energies, and exponential decay 
into the lower energies (which will again be mollified from photons descending from higher 
redshifts). However, for this case the position of the kink will be fixed at around 100 MeV, 
independent of the reheat temperature. 

Let us compare the above prediction with the recent COMPTEL and EGRET data. The 
EGRET collaboration found that the flux in the energy range 30 MeV — 100 GeV is well fit 



by the single power law [|28[ 



dUr 



- ^— (cm 2 sec sr MeV)" 1 , (34) 

451 MeV/ v ; y J 



while the COMPTEL data |2jJ in the range 0.8-30 MeV can be fit |3 to the power law 

d J ( uj \ ~ 2 - 38 

6.40 x 10~ 3 - (cm 2 sec sr MeV)" 1 . (35) 



dcu VI MeV, 

Below 0.8 MeV there is large increase in the measured flux. Thus, the best bounds are 
found by comparing the measured flux to predicted flux at Wkink or at 0.8 MeV, whichever 
is larger. Because of the rapid rise of the predicted spectrum relative to the measured 
spectrum, a change in the kink position can change the bound on n on the order of 0.01, 
which we consider within the accuracy of our calculation. The bounds on n from the diffuse 
gamma-rays are specified by the dot-dashed line in Fig. |]. The bound terminates when all 
but the exponential tail of the PBHs decay prior to a redshift of 700, since the optical depth 
at such early times exceeds unity, as discussed above. 

We may compare our results to those derived by Yokoyama [|31||] where the author 
placed bounds on mass fraction of PBHs at /3(tR#) = fi(£##), using the initial mass 
function, Eq. (|i"0p. He found that the bounds on /3 did not differ significantly from the 



4 After this work was completed, we became aware of Ref. |32j that also utilized the critical 
collapse initial mass function to derive bounds on the PBH mass density by requiring that LSPs 
(in supersymmetric models) not be overproduced. 
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previous bounds derived using the standard initial mass functions, except for the bounds 
coming from diffuse gamma-rays. In the latter case, applicable for horizon masses in the 
range Mh > 5 x 10 14 g, he found more stringent constraints. Our bounds on n, translated 
into bounds on (3, agree with his bounds coming from energy density constraints except for 
the case of larger reheat temperature, since we included the proper scaling of the energy 
density of photons emitted after PBH decay. Thus our bounds on (3 can differ by many 
orders of magnitude. Our bounds coming from diffuse gamma-rays can also differ by orders 
of magnitude, but in this case it is for a different reason. Yokoyama determined his bound 
on (3 by imposing the constraint on Q p bh{to) derived in Ref. 0. However, when we change 
the initial mass function we also change the diffuse gamma-ray spectrum significantly in 
both shape and normalization, as discussed above. Thus, it is inappropriate to directly take 
the bounds from Ref. and apply them to the case with the new initial mass function, 
Eq. ([TOD. We find that our bounds on f3 from diffuse gamma-rays are more stringent than 
those determined in Ref. |HJ in the range Mh > 5 x 10 15 g by several orders of magnitude. 



V. ROBUSTNESS OF THE BOUNDS 

Let us consider the robustness of the bounds. We might worry that the bounds are 
highly sensitive to the choice of parameters given the sharpness of the initial mass function. 
Indeed, in the case where it is assumed that the PBH mass is given by Eq. (0), the bounds 
are n are exceptionally sensitive to the exactness of this relation. This is clear from the 
exponential factor in Eq. ([32|). Given the initial mass function calculated by Jedamzik 
and Niemeyer, Eq. (|i~0|), we must check the sensitivity to the parameters S c , k and ctq. In 
Ref. fL6f , the authors tested the scaling relation (|3]) using several different initial shapes 



density perturbations shapes. They found (6 C = 0.70, k = 11.9), (5 C = 0.67, k = 2.85), 
(5 C = 0.71, k = 2.39), for Gaussian, Mexican Hat and fourth order polynomial fluctuations, 
respectively. 

We varied the value of 5 C between 0.60 — 0.80 and found that the bounds changed by 
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at most 0.01. The sensitivity to the variation being maximal at the smaller values of Trh. 
Given that the initial mass function is peaked at a number smaller than kM^, the sensitivity 
is increased at smaller Trh because a is an decreasing function of Trh- Variations in k are 
equivalent to a scaling in Trh- Thus, varying k by an order of magnitude has essentially no 
effect on the bound. Lastly, let us consider the sensitivity to the parameter a . The value 
we used for uq in Eq. @ was calculated in Ref. flTj using the result [|HH F| 



5 = 1.91 x 10 - 5 exp[1.01(l-n)] 

where r is a measure of the size of the tensor perturbations. The la observational error being 
7%. The fit, Eq. (p6[) , is good to within 1.5% everywhere within the region 0.7 < n < 1.3 



and < r < 2. The authors of Ref. |19| quote a 9% uncertainty in Eq. ( p6|) at la, once 
uncertainties in the systematics and variations in the cosmological parameters are taken 
into account. The value in Eq. (§) was determined ignoring tensor perturbations. Given 
that o"o scales with 5q we find that varying a® at the 2a level has no effect on our bound at 
the level of 0.01. On the other hand, including some contribution from tensor perturbation 
will weaken the bound. We found that taking r = 2 weakened the bound by 0.01 — 0.02 
throughout the range in the reheat temperature. 

We can also consider the effects of a non- vanishing Q^. Bunn et al. [19] extended their 
results to this case and found 



r I 1 ni -.s m-5 exp[1.01(l n) n -0.80-0.051ogf7o /-\ , nis/ t\ n no n ^ 

oo|fi A = 1-"1 x 10 Uq s (1 + 0.18 [n — 1) — O.Odr^Aj • 

+ (0.75- 0.13 tt 2 A )r 

(37) 

If < r < 2, we can express 8o\n A extracted assuming a nonzero cosmological constant to a 
very good approximation by a scaling of 5q extracted without a cosmological constant 



5 It should be emphasized that this result assumed the spectra can be approximated as a power 
law over the range of k that COBE probes. We are then making the further assumption that n is 
constant down to the mass scales of relevance for PBHs. 

15 



8 \u A - Q ^-^^ 6o , (38) 

(where Q + = 1) and thus <Jq also acquires a correction. Consequently, the bound on n 
is shifted by 

2 (-0.8 -0.05 In fi ) lnfi , Qn , 
A "^-"I"a= 42.9 + ln(TWl0 8 GeV) " (39) 

In Fig. ||| we show the above correction as a function of Trh for several choices of If we 

take ~ 0.7 as recent observational data suggests, our bounds on n strengthen by about 

0.03 — 0.06 for Trh between 10 16 — 10 3 GeV respectively, as shown in Fig. f|. 

We can also calculate bounds on the mass variance at reheating |33[] which is essentially 

model-independent. If the relation Eq. (|7|) is violated by, for example, a power spectrum 

with a spectral index that depends on scale, then our previous bounds on n cannot be 

applied. However, given a inflationary model one could in principle calculate the power 

spectrum, normalize to the COBE data at our present epoch, and then match onto the mass 

variance at reheating. In Fig. |^ we show the bounds on the mass variance from both the 

density bounds as well as the bounds from the diffuse gamma-ray observations. Notice that 

the diffuse gamma-ray observation bounds on cx(M#) are a significant improvement over the 

density bounds in the applicable range of reheating temperatures. 



Finally, we must address the issue of non-Gaussianity. It has been pointed out f34 
that skewness could very well be important for PBH formation given that its effects are 
amplified in the tail of the distribution P[5], which contributes to PBH formation. In 
general, the amount of non-Gaussianity expected is highly model dependent. Bullock and 
Primack investigated several inflationary models to study the amount of non-Gaussianity 
one would expect at larger values of 5. They calculated P[5] for three toy models, and 
found in one case no deviation from Gaussianity and in the other two found a significant 
suppression in the probability of of large perturbations. However, as was pointed out in 
Ref. [|1(|, while these effects can drastically effect the PBH mass fraction (3, we expect that, 



even in the most extreme case considered in Ref. , the effect on the bound on n is only at 
the level of 0.05. For hybrid inflation, where the approximation that n is constant actually 
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holds, the perturbations are in fact Gaussian due to the linear dynamics of the inflaton field 
f35|j . Therefore, these bounds should be applied to specific models, with the roughness of 
the bound determined by the deviations away from Gaussianity. 

VI. CONCLUSIONS 

We have calculated the density of primordial black holes using the the near critical 
collapse mass function that results in a spectrum of PBH masses for a given horizon mass. 
The normalization of the PBH mass spectrum was determined using the COBE anisotropy 
data that allowed us to set bounds on the spectral index n as a function on the reheat 
temperature. We find that restricting the density of PBHs to be less than the critical 
density corresponds to the restriction that the spectral index n be less than about 1.45 to 
1.2, throughout the range of reheating temperatures resulting after inflation, 10 3 to 10 16 
GeV respectively. (The precise limits are shown in Fig. [I].) For a smaller range of reheating 
temperatures, between about 10 7 to 10 10 GeV, significant PBH evaporation occurs when the 
optical depth of the universe is less than one. Hence, we found a slightly stronger bound on 
the spectral index by restricting the cosmological PBH evaporation into photons to be less 
than the present-day observed diffuse gamma-ray flux. Due to the extreme sensitivity of the 
PBH mass density to the spectral index, effects such as the indirect photon flux from PBH 
evaporation into quarks and gluons which fragment into pions or the formation of a QCD 
photosphere are completely negligible when calculating the bound on n. We should also 
remark that slightly stronger bounds on n for larger reheating temperatures > 10 10 GeV are 
expected from PBHs that decay during the epoch of nucleosynthesis. 

If the universe is vacuum-energy dominated, there are corrections to our bounds on n 
that can be substantial. We calculated these corrections for a range of Q\ and applied them 
to our bounds on n for the case of Q\ = 0.7. The improvement is apparent by contrasting 
Fig. |T| with Fig. Finally, we calculated bounds on the mass variance at reheating. These 
bounds in principle could be used to constrain any given inflationary model, once the power 
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spectrum is calculated. 
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FIG. 1. The upper bound n as a function of Trh by requiring that Q p bh,evap(to) + ^pbh(to) < 1 
(solid line), & p bh(tRH) < 1 (dotted line), and that the PBH photon spectrum does not exceed the 
diffuse gamma-ray background (dot-dashed line). If PBHs leave a Planck mass relic, an additional 
bound is present for large Trh (dashed line). Note that the endpoint of each line within the figure 
is where we elected to stop calculating the bound on n, due to the presence of another stronger 
constraint as shown. 
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FIG. 2. Examples of the direct cosmological photon spectrum from evaporating PBHs with 
n « (1.340,1.309,1.285), T RH = (10 7 , 10 8 , 10 9 ) GeV for the (top, middle, bottom) solid lines. 
Notice that the spectrum scales roughly as dJ/dujQ oc o;^ 3 for ujq > uJkink- At low energies 
(ujq <C 100 MeV), the spectrum is modified by the production of quarks and gluons emitted from the 
PBHs that fragment into pions, which then decay into photons. At high energies (loq > 100 MeV), 
the spectrum is modified if a photosphere forms around the PBH. Since the normalization of the 
PBH mass spectrum is extremely sensitive to the spectral index, both these effects are completely 
negligible when calculating the bound on n. 
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Trh [GeV] 

FIG. 3. The correction An = n — n\n A as a function of Tjm to all our bounds in Fig. |l] induced 
by assuming a cosmological constant J7a = (0-9, 0.7, 0.5, 0.3, 0.1) for the five lines in the figure from 
top to bottom respectively. 
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FIG. 4. The upper bound on n as a function of Trh as in Fig. |l|, except that J7a = 0.7. 
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FIG. 5. Same as Fig. ||, except that we show the upper bound on cf{Trh) as a function of Trh- 
This bound is independent of the spectral index n. 



24 



